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Abstract 

In this paper we suggest an approach to analyse 
the motion of a test particle in the spacetime of 
a global monopole within a /(i?)-like modified 
gravity. The field equations are written in 
a more simplified form in terms of F(R) = 
■ Since we are dealing with a spherically 
symmetric problem, F(R) is expressed 
radial function T[r) = F(R(r)). So, the choice 
of a specific form for f(R) will be equivalent 
to adopt an Ansatz for J-(r). By choosing 
an explicit functional form for J~(r) we obtain 
the weak field solutions for the metric tensor, 
compute the time-like geodesies and analyse 
the motion of a massive test particle. An 
interesting feature is an emerging attractive 
force exerted by the monopole on the particle. 
PACS numbers: 04.50. + h, 04.20. - q, 
14.80.Hv 



1 Introduction 

A global monopole is a kind of topological 
defect which arises in certain gauge theories 
due to the spontaneous symmetry breaking 
(SSB). Such SSB processes may have been 
undergone by our universe as a consequence 
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of many phase transitions in the early times. 
Global monopoles arise in a theory exhibiting 
SSB of the global gauge group 0(3) to U(l) 
[U [2]. The gravitational field of this object 
in the context of General Relativy (GR) was 
investigated by M. Barriola and A. Vilenkin 
[3]. There, the authors have shown that 
the spacetime associated with this object 
is characterized by a non-trivial topology 
represented by a deficit solid angle and presents 
a non- vanishing scalar curvature. Later, Barros 
and Romero [3] analysed the gravitational field 
of a global monopole in the context of the 
Brans-Dicke theory of gravity considering the 
weak field approximation. There a comparison 
is made with the corresponding result obtained 
from GR. 

In the late ten years, several types 
of modified theories of gravity have been 
suggested as possible alternatives to explain 
the late time cosmic speed-up experienced by 
our Universe, one of these theories are the so- 
called f(R) theories of gravity [5j EJ [7]. Such 
theories avoid the Ostrogradski's instability 
which is commonly observed in general higher 
derivatives theories [H [9] . 

The gravitational field of a global monopole 
in such modified theory of gravity has been 
investigated in [10]. There, we have found 
solutions for the metric tensor in the weak field 
approximation considering a specific Ansatz for 
the functional f(R)- The field equations are 
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expressed in terms of the function F(R) = 
d ^j$ , in a similar approach as developed in [Tl] 
and |12j . Since the system under investigation 
has a spherical symmetry, it was possible 
to write F(R) as a function of the radial 
coordinate, r, only. In the present paper we 
intend to follow the same procedure previously 
adopted by us, nevertheless going further and 
investigating a wider class of arbitrary n-degree 
polynomial functions J-(r). We consider some 
approximations. For instance, the weak field 
limit and the modified gravity as a small 
correction on GR. 

In the scenario described above, we analyse 
the motion of a massive test particle in the 
spacetime of a global monopole from a classical 
point of view. The analysis of the classical 
motion of a massive test charged particle 
in the spacetime of a global monopole has 
been developed in [13J, taking into accout 
the presence of the induced electrostatic self- 
interaction. A study of the motion of particles 
in the field of a non-Abelian monopole in 
Einstein gravity was made in |14j . 

This paper is organized as follows: In 
the section [2] we review briefly the model 
proposed by Barriola and Vilenkin for the 
global monopole in the context of GR. In the 
section [3] we introduce the modified theory of 
gravity and the field equations in the metric 
formalism, which is a crucial point to obtain in 
section U] the solution for the global monopole 
in f(R) gravity. In the section [5] we analyse 
the classical motion of a test particle in this 
spacetime. Finally, we leave for section [6] the 
most important concluding remarks about this 
paper. 

2 A Brief Review about 
Global Monopole 

The model describing a global monopole is 
given by the following Lagrangian density [3]: 

c = \d^d»4> a -\\(r<t> a -v 2 ) 2 • (i) 

Where <j) a is a self-coupling triplet of scalar 
fields and a = 1,2,3. Plugging the 
energy-momentum tensor arising from (pTJ) into 



the Einstein equations, spherically symmetric 
solutions can be obtained, by adopting for the 
line element and the matter fields the Ansdtze 
below: 

ds 2 = B{r)dt 2 - A{r)dr 2 - r 2 d£l 2 (2) 

and 

r = vKr)y, (3) 

with x a x a = r 2 and dti 2 = d6 2 + sin 2 9dip 2 . 
Barriola and Vilenkin consider points far 
from the monopole's core where the energy- 
momentum tensor associated with the static 
matter field © can be approximated as 

Z?«diag(|£ > £ > 0,0) . (4) 

Taking these assumptions into account they 
have obtained the following solutions for the 
functions B(r) and A(r) 

B = A" 1 « 1 - 8ttGt] 2 - 2GM/r , (5) 

being M the mas of the monopole. 

Furthermore, they have computed the light 
deflection performed by the gravitational field 
of the monopole and found 

50 = 8TrGri 2 l(d + l)- 1 , (6) 

where d and I are the distances from the 
monopole to the observer and to the source, 
respectively. The light ray is constrained to lie 
at the equatorial plane 6 = \. 

3 The f(R) Gravity in the 
Metric Formalism 

In a f(R) theory, the action associated with the 
gravitational field coupled to a matter field is 
given by: 

S = h I dix ^9f{R) + S m , (7) 

where f(R) is an analytical function of the 
Ricci scalar, R, k = 8ttG, being G the Newton 
constant, and S m corresponds to the action 
associated with the matter fields. By using the 
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metric formalism, the field equations turn out 
to be: 

1 



F{R)R^ - -f(R)g^-V^uF(R) 

+ g^nF(R) = KT™, (8) 

where F(R) = and T™ is the standard 

minimally coupled energy-momentum tensor 
derived from the matter action, S m . 

By taking the trace of equations ([8]) we 
get another one which shows us explicitly 
an emerging scalar degree of freedom in the 
modified gravity: 



F{R)R - 2/(12) + 3DF(R) = kT t 



(9) 



This expression allows us to write the function 
f(R) in terms of its derivatives as follows 

f(R) = i (F(R)R + 3DF(R) - kT 171 ) . (10) 



When (|T0|) is plugged into ([8]) it is possible to 
express the field equations in terms of F(R) as 
shown below: 

F(R)R^ U - V^X7„F(R) - kT™ = 

^ [F{R)R - DF(R) - KT m ] . (11) 

Of course, it is simpler to handle diferential 
equations involving F(R) than to deal with 
them in terms of f(R), since the former 
function is the first derivative of the latter 
one. So when this replacement is made the 
order of the field equations gets reduced by 
an unity. Other important reason for working 
out our approach looking always at F(R) as 
a dynamical variable, instead of f(R), is the 
new feature arising from equation (jlOf) . This 
equation says that a further scalar degree of 
freedom appears in in a /(12) theory and, 
according to that equation, is exactly F(R) 
which carries this information. 

From the expression fill j) we can see that the 
combination below 



C u 



F(R)R^-V tl V ll F(R)-KT^ l 

9pLfl 



, (12) 



with fixed indices, is independent of the 
corresponding index. So, the following relation 



Cu — C u 







(13) 



4 The Field Equations 

Substituting the energy-momentum tensor 
(J2]) into (fT3j) . in [10] we have derived 
the field equations associated with a global 
monopole system in modified theories of 
gravity. Considering for the metric tensor the 
Ansatz given in (J2]) we have: 



, ,B' A'\ ( R' A' 



= 0, 
(14) 



and 



p p 
4B + 4AB - ArB— + 2r 2 B'— 
F T 

r ^B' + 2r 2 B 



B A, 
(B' A!\ AABkt? 
+ 2Br [-B + A) T^ = °' 

(15) 

where we have expressed F(R) as F(R(r)) = 
F(r), since the problem we are analyzing has a 
radial symmetry, and primes mean derivative 
with respect to r. Those equations may be 
simplified if we consider the following definition 



a B ' A ' 

~B + ~A 



(16) 



which will allow us to write the field equations 
as follows 



/3 _F" IF' 

7 ~ T ~ Yf ' 



(17) 



and 



77/ 77/ 
AB + AAB - ArBy + 2r 2 J B / — + 2r 2 B" 



r 2 B'p + 2BrP- 4A f r > =0. 



holds for all \x and v. 



(18) 



As it has pointed out in [TQ, it is possible to 
verify that if one plugs cu = (Brans-Dicke 
parameter) into the equation (9) of reference 
[1] and replaces 4>{r) by F(r), the equation 
(]17p is recovered. Which is a clear evidence 
of the equivalence of the f(R) gravity in the 
metric formalism with a Brans-Dicke gravity 
possessing a parameter oj = 0, as previously 
discussed in |18| . Furthermore this fact offers 
us a deep ground to work with the field 
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equations in terms of F(R), since the results 
obtained by us can be compareable with those 
ones found via Brans-Dicke gravity. In this 
case, the role of the scalar degree of freedom in 
f(R) gravity is played by the function F(R). 

In order to simplify our analysis possible, we 
shall assume some approximations. First of 
all, we choose to handle the field equations in 
the weak field regime, which means to consider 
that the components of the metric tensor are 
given by: B{r) = 1 + b(r) and A(r) = 1 + a(r) 
with \b(r)\ and \a(r)\ much smaller than unity. 
We also consider that the modified gravity is 
just slightly deviation from GR, i.e F{r) = 
l + tp(r), with | ip(r)\ « 1. Furthermore, since 
Grj 2 ~ 10 -5 in G.U.T., we may keep only terms 
linear in Grj 2 and ip(r). By ruling out all the 
crossed terms of such small quantities and their 
derivatives, we get a linearized form for the 
original field equations 



and 



r 



(19) 



and 



2a - 2rtp' + r(3 + r 2 b" 
- 2kt] 2 = . 



(20) 



Closed results for equations (|T9[) and (|20|) 
can be obtained by adopting for F(r) a specific 
Ansatz which means merely to solve them for 
a given ip(r) as it will be done in the next 
subsection. 

4.1 The Solution for F(r) = 1 + ip r n 

Let us consider for the function ip{r) a power 
law-like Ansatz, namely ip(r) = V'o 7 ''™) where 
ijjQ is a constant parameter to be determined 
associated with the modification of the gravity. 
It is convenient to impose ip(r) as being regular 
at the origin which implies that ift(Q) = 0. 
This assumption automatically rules out all the 
negative powers of n. Moreover, we assume 
that ip{r) is an analytical function of r, so 
that it admits a Taylor series expansion. The 
equations (|19p and (|20f) for the aforementioned 
Ansatz read, respectively 



a + b' = n(n — l)tpo rr 



2a - 2nip r n + r(a' + b') + r 2 b' 
-2kt) 2 = . 



(22) 



In order to ensure that the metric is 
asymptotically flat as ipo — > 0, the integration 
constant arising from pip is set to be zero, 
therefore a following relation between a{r) and 
b(r) is obtained: 



a + b = (n — l)ipor n . 



(23) 



The relation above can be used in (|22p in order 
to express it in terms of b(r) 

r 2 b" - 2b+ ( n +l)(n- 2)ip r n -2kt] 2 = , (24) 



whose solution is 

b(r) = — — C2r 2 — ipo 1 "™ ~ K V 2 
r 



(25) 



where the integration constants c\ and C2 are 
defined, respectively, as c% = —2GM and 
C2 = 0, to recover suitably the Newtonian 
potential present and due to the absence of 
a cosmological constant in the model we are 
analysing. Since B{r) = 1 + b(r) we have: 

B{ r ) = 1 - — 8ttGi] 2 - Vor" . (26) 



(27) 



From equation (|16p we obtain 

A(r)B(r) = aoe^' 1 ^^ 



where we set the integration constant ao to be 
unity in order to have a spacetime obeying the 
asymptotical flatness condition. Therefore, 



A(r) 



3 (n— l^or" 



2GM 



87rGy - Vor r 
(28) 

Following the same reasoning of Barriola and 
Vilenkin, we drop out the mass term which 
is negligible at astrophysical scale, so we may 
write 

B{r) = 1 - MGif - i; r n (29) 



and 



(21) 



A(r) = e ( n - 1 ^ or " (1 - 8vrCV - ^ r n ) 1 



(30) 
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A binomial expansion could be applied in the 
solution ([30]) so that it would be expressed as 



A(r) wl + MGrj 2 + nip r r 



(31) 



It is important to recall that our analysis 
is restricted to a specific range of the radial 
coordinate r, which is 



5 < r < 



1 



(32) 



, -l 



being 5 ~ (A?] 1 / 2 ) * of the order of magnitude 
of the monopole's core. The reason for that, 
is because we have derived the field equations 
outside the monopole's core where h ~ 1; 
moreover, due to the weak field approximation 
we must have \ipo r n \ < 1. 

The solutions ([29]) and ([3T]) above allow us 
to compute the corresponding Ricci scalar 

R = - [( n - l)( n - 2) + 2(n + 2)] ^ n " 2 

- ^ ■ (33) 

from which it is possible to determine the 
explicit form of f(R). The procedure to 
be followed consists in inverting the equation 
above by expressing r as a function of R, 
then plugging r = r(R) into J-{r) = 1 + 
%[)Qr n and finally performing the integration 
of F(R) which will give f(R) plus an 
integration constant, which we may discard if 
no cosmological constant is taken into account 
in the theory. For any physically viable 
f(R) theory under consideration, the following 
stability conditions must be fullfilled [T5]-[T8]: 

• d l^P > (no tachyons); 

• d fjn' > (no ghosts); 



lim R ^oo tj = and lim^ 
is recovered at early times) 



dA 

00 dR 



(GR 



where A = A(R) is defined as A = f(R) - R. 



GR. Let us consider the transformation of 
coordinates below: 

B(r) = p(r) (l - 8irGr] 2 ) , (34) 
A(r)dr 2 = p(f) (l + 8irGr] 2 ) df 2 , (35) 
r = p l l 2 (f)r , (36) 

where p(f) is an arbitrary function of f to be 
determined and p{f) = l + q(f) with \q(f)\ < 1. 
Diferentiating equation (|36[) we have 



dr 2 



dq 
dr 



df' 



(37) 



Substituting Eq. (I37jl into (135]) we obtain the 
following result for q(r) (we keep only linear 
terms in q(r), ipor and Grj 2 ): 



then 



p{r) = 1 — ipor n . 



Thus we can write the line element 
coordinate r as follows 



(38) 

(39) 
in the 



ds 2 = (1 - i> f n ) [(1 - 8^Gt? 2 ) dt 2 
- (1 + 8ttGti 2 ) dr 2 



~" 2 [dO 2 + W6d(t) 2 )] . 



(40) 



If we rescale the time coordinate and redefine 
the radial coordinate as r = (l + AttGt] 2 ) f we 
arrive at the line element below 



ds 2 



(1 - ip r n ) [dt 2 - dr 2 - (1 - 8tt<V) 



xr 



+ sin 



edcf 2 )] . 



(41) 



An important feature arises if we are 
interested in analyzing the deflection of light 
in this metric. As it is well known, the 
deflection angles are always preserved for two 
metrics related by a conformal transformation. 
Therefore, the deflection of a light ray by 
the monopole in the present modified gravity, 
considering weak field approximation, will be 
the same of that one previously obtained in 
Ref.[3]. 



4.1.1 The conformal relation with 
Barriola-Vilenkin metric 

It can be verified that the line element 
described by ([29]) and ([3T]) are conformally 
related to the global monopole solution in 



5 The Classical Motion of a 
Test Particle 

In this section we analyse the classical motion 
of a test particle in the spacetime whose metric 
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tensor has the form ([5]) with components given The first and second conditions provide us the 
by (|29p and ()31 j) . From the corresponding line following system of equations: 
element we can define the Lagrangean of a test ^_, 2 j2 

~B~~r 



particle moving on this geometry as follows: — - — 1 = (47) 



-V = 2C = B(r)i 2 -A{r)r 2 and 
drj q 

- r 2 9 2 - r 2 sin 2 6ip 2 , (42) ~q t 



B r2 1 



. (48) 



where the dot means derivative with respect to Whose solution for E and L are, respectively, 
t. given by: 

For orbits at equatorial plane, i.e., with _ / 2 

# = the corresponding canonical momenta, V 2B — B'r 

Pa = j^k are, and 



B'r 3 

p t = E = B(r)i, Pe = -L e = 0, L = \J w _ B , r • (50) 

Pr = — A(r)r , = —L<p = —r 2 (p ,(43) The angular velocity for any equatorial orbit 

, , - . _ , , is defined as 

where the constants of motion E and L v 

are interpreted, respectively, as the energy, ^ = ^_ _ dtp/dr 

and angular momentum, per unit mass, in ip- dt dt/dr 

direction of the particle. SO; from the equations ggj and we obtain: 

For a massive particle the relation below 
holds B L [B 7 

dx»dx v 1 , tA . Q = ^v = \hr- 52 

9^ W = 1, (44) r*E V2r 

t i - , i , ■ It is well known that the velocity of the test 

from which we obtain -L 

particle along the i-direction is [19, 20J: 

[l + (n-l)^ r n }r 2 + V eff (r)=E 2 . (45) < { 

v 2 = _^dx^_ (53) 

That shows a position-dependence of the 500 dt 

particle's mass. Moreover, V eff (r) is the ThuS) the tangential velocity v v reads 

effective potential energy associated with the 



test particle and is defined as ^2 _ 9w (dip\ _ g w , o2 



v z = rr = __f!£o 2 . (54) 

v goo \dtj g o 



r 2 ) ' By working out the expression above and 



V eff (r)= {1 - SttGt] 2 - 4> r n ) ( 

(46) keeping only linear terms in Grf and ipo we 

with L = Lp. obtain 

In the next subsections, we shall analyze v _ , /_i n ^, or n ^55^ 

more detailedly the classical motion of the test » 2 

particle. The expression above says that circular orbits 

will be physically allowed only if ipo < 0. 

5 1 The tangential motion Hence, such equation can be rewritten as 

For this case it is necessary to require the stable /I 1 „ /r - c \ 

. , .. . , = V n n lw r n , (56) 

circular motion conditions given by \ 2 

which shows us that in the case under 
consideration, with the mass of the monopole's 



(i) f = , 



(ii) Q 5 r 5 = , core considered as negligible, the circular 

motion of the test particle is consequence of 

d 2 V tt 

(iii) — > . the modification parameter of the gravity, Vo- 
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5.2 The emerging extra force 

One important physical property observed 
when the global monopole is analyzed within a 
f(R) gravity is the gravitational force exerted 
by the monopole on a massive particle moving 
nearby. It is also important to emphasize 
that this feature is absent in GR, since the 
Barriola-Vilenkin monopole has a metric with 
<7oo =cons t- which means that no gravitational 
force exist. On the other hand, based on our 
previous result, we have found 



(57) 



which gives rise to a radial force acting on the 
particle for ipQ ^ 0. 

As it is well known, the motion of the 
test particle in a weak gravitational field is 
described by the equation [T9] 



ldh 



x 



00 



2 dx l 



(58) 



being Hqq the deviation form the unity in 
<?oo- In order to obtain the force given by 
(f58|) . let us express the metric ([2]) in galilean 
coordinates, g^ v = rj^ + h^, by performing 
the transformation below 



t = (1 - AttGt] 2 ) T , 



(59) 



and 



,-ii-4 rt y-* , °^ h «*'*'V, 



n 



(60) 



which will imply 



ds z 



(1 - 8irGr/ 2 - ip K n ) dT 2 
[1 - WirGrj 2 - ^j TZ n 
87rGr] 2 ln(|V> | ft") ' 
n 

(dK 2 + K 2 dn 2 ) , 



where 1Z = \J x 2 + y 2 + z 2 . Hence, by working 
out (EH) for (EH) we obtain 



ntpoTZ 71 - 1 x i 



whose corresponding force experienced by an 
unit mass particle is 



n-l 



-K . 



(63) 



Because of the negativeness required for ipo 
such force will have an attractive nature which 
may be evidenced by expressing (f63l) as 

-+ n\^\TZ n - 1 



-K . 



(64) 



Therefore, as a consequence of modifying GR 
we have obtained an extra and attractive force 
exerted by the global monopole on a massive 
test particle. 

It is useful to investigate what are the 
possible motions that the particle can perform. 
In order to do that, let us express the effective 
potential in terms of a set of dimensionless 
variables defined as follows 



( V X 1/2 )r, |^o| 



L~2= L 2 ( V 2 X) 



a 



|*>| (riX 1/2 
= 1 - 8ttGt] 2 . (65) 



With these changes the range of validity ([32 
is redefined as 

1 



1 < f < 



l^o I 



(66) 



A remarkable feature which is important to 
be emphasized is that the range of validity 
presented above is shorter the higher is the 
value of n. We mean, the size of the region in 
space in which some important change will be 
observable due to the modification of gravity 
tends to zero as higher powers of ip(r) are 
considered. 

In terms of those new variables the effective 
potential (|46p may be written as 



n 



V eff (r) = a 2 + \j \r n + a 2 ^ - ^L 2 r n ' 2 . 

(67) 

By assigning numerical values for the 
parameters a 2 , ipo, L and n, present in the 

(61) 

equation above, we can analyse the behavior of 
the possible profiles of the effective potential. 
Let us make plots of V e ff(f) against f for 
different values of tpo and then observe how the 
effective potential behaves for especific values 
(62) of n, namely n = 1, 2 and 3, respectively. 
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5.2.1 Plot for n = 1 

A plot of V e ff(f) against f is exhibited in the 
Fig 1. We have sketched the graph by assigning 
specific values for a 2 and ipo. In this plot two 
distinct configurations of effective potentials 
are presented: In the one represented by solid 
line, compatible with the context of modified 
gravity, we can see that the attractive force 
traps a test particle moving with energy greater 
than the minimum potential. In the second 
one, represented by dashed line, the parameter 
responsible for modifying the gravity ipo is 
set equal to zero and the potential has no 
extremal. So, in the latter, the particle will 
not be trapped by the monopole, as predicted 
by Barriola-Vilenkin since this monopole does 
not exert any gravitational force on the massive 
particles moving nearby. 




10 20 30 40 



Figure 1: Plot of the effective potential V e ff vs f for 
two different ,set of values of a, ipo and L. The solid line 
shows a minimum of the effective potential for a 2 = 0.9, 
ipo ~ —0.02 and L = 1, which means that a particle 
moving with energy E 2 > V e //(r min ) will be trapped 
by the monopole. On the other hand, the dashed line 
expresses the corresponding effective potential within 
GR, with a 2 = 0.9, ipo = and L = 1. In the latter 
case, the potential has no minimum and the particle 
cannot be trapped by the monopole as expected from 
the result of Barriola-Vilenkin. 

5.2.2 Further plots for n = 2 and 3 

The next examples we intend to analyse are 
that ones concerning to n = 2 and 3. The plots 



V e fj vs r shown below represent the profiles 
of the effective potential (|67|) for each one of 
those values of n. The numerical values of ipo 
adopted in each case were ipo = —0.02. The 
different scales of the horizontal axis in each 
plot remarks the dependence of the range of 
validity ([661) upon the value adopted for n. 




2 3 4 5 6 7 8 



Figure 2: Case n = 2. For 
such case the range of validity is 
1 < f < 7.07, so it is reasonable 
to leave the radial coordinate to 
vary between 1 and 8 as it is well 
expressed in the plot. 




12 3 4 



Figure 3: Case n — 3. Now 
the range of validity is 1 < 
f < 3.68, which suggests the 
variation shown in the horizontal 
axis, namely between 1 and 4. 

As we can see, for increasing values of 
n, the graphs trend to keep similar profiles, 
always having a well-defined minimum which 
is necessary for the arising of the extra force. 
Since the force (j63|) is attractive, it will be 
possible to observe a trapping of the test 
particle by the Global Monopole if the particle 
moves with energy greater than the minimum 
potential. In each case sketched above there 
will be two turning points for the motion of 
the particle. This trapping is a feature absent 
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within GR, being an exclusive consequence of 
this modification in gravity we did. 

6 Conclusion 

In this paper we have analyzed the classical 
motion of a massive test particle in the 
gravitational field of a global monopole in 
the f(R) gravity scenario in the metric 
formalism. By this formalism, this modified 
gravity contains a massive scalar degree of 
freedom in addition to the familiar massless 
graviton and it turns out to be equivalent to 
a Brans-Dicke theory. In order to simplify 
our analysis we have considered solutions in 
the weak field approximation, which implied 
that we are considering a theory that is as a 
small correction on GR. The latter condition 
was explicitly considered by assuming J-(r) = 
1 + ipor n , for IV'o 7 ""! < 1) being n an positive 
integer number. 

Following the above mentioned approach, 
the solutions found by us correspond to 
small corrections on goo an d gu components 
of the metric tensor, only. Being these 
new components given in (|26|) and (|28D . 
respectively. From the results obtained, we 
can observe that, a small correction on the 
Ricci scalar, given now by (|33|) . takes place. 
Moreover, we have also verified that these 
solutions are conformally related to that one 
previously obtained by M. Barriola and A. 
Vilenkin as shown in (j41|) , what ensures us that 
the deflection of light in these two spacetimes 
will be the same. 

As to the classical motion, we have derived 
the differential equation, given by Eq. (|45j) . 
which allows us to analyse separately the 
tangential and radial motions. As our main 
conclusion, we have observed that the presence 
of the parameter ipo is essential to provide 
stable circular orbits for the particle. This 
is a new feature, because in the context of 
general relativity and discarding the mass term 
in the metric tensor, there is no Newtoninan 
gravitational potential produced by a global 
monopole, consequently, no stable orbits for 
the particle is possible. 

An important result we have found was 
an emerging force arising within the modified 



gravity, which depends on the parameter ipQ 
and also on the values attributed to n. Such 
force, as verified by us, has an attractive nature 
and becomes stronger for larger values of n. 
This is a new effect arising in the modified 
theory, being absent within Barriola- Vilenkin 
model. In the last section we have sketched 
graphs for the effective potential for specific 
values of n, namely n = 1, 2 and 3 and 
observed that their behavior are similar, in 
the sense that they always present attractive 
emerging forces. 
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